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Neutrino oscillations are described and interpreted using the front form of Hamiltonian dynamics 
in the Feynman-Gell-Mann-Levy version of an effective quantum field theory in which leptons 
interact with whole nucleons and pions instead of quarks. The interactions are treated in the 
lowest-order perturbative expansion in the coupling constants Gf and F„ in the effective theory, 
including a perturbative solution of the coupled constraint equations. Despite missing quarks and 
their binding mechanism, the effective Hamiltonian description is sufficiently precise for showing that 
the standard oscillation formula results from the interference of amplitudes with different neutrinos 
in virtual intermediate states. This holds provided that the inherent experimental uncertainties 
of preparing beams of incoming and measuring rates of production of outgoing particles are large 
enough for all of the different neutrino intermediate states to contribute as alternative virtual paths 
through which the long-base scattering process can manifest itself. The front-form result is the same 
as the previously derived instant-form result despite that the long-base scattering process is traced 
in space-time differently and the relevant interaction Hamiltonians are constructed differently in the 
different forms of dynamics. 



I. INTRODUCTION 

Contemporary theory of neutrino oscillations can be 
developed using different forms of relativistic Hamilto- 
nian dynamics. A priori possible forms had been classi- 
fied by Dirac over 60 years ago He generally distin- 
guished three forms. In the first of these forms, which 
is used most commonly and which Dirac called the in- 
stant form (IF), a Hamiltonian H generates the evolu- 
tion of a system in time t of some inertial observer. This 
means that the evolution is traced in terms of data that 
change from one space-time hyperplane of constant t to 
another. These hyperplanes are called instants. In the 
second form, called by Dirac the point form (PF), the 
evolution is traced from one hyperboloid in space-time 
to another. The PF distinguishes a point in space-time 
and the operators of spatial momentum in it, or gener- 
ators of translations in space, involve interactions. This 
is why the PF is not popular despite that the Lorentz 
symmetry is represented in it purely kinematically. In 
the third form, called by Dirac the front form (FF), the 
evolution of a system is traced from one space-time hy- 
perplane of constant x + to another, where one typically 



uses the convention x + 



x along with the nota- 



tion ± = ± 3 for indices of all tensors. A hyperplane 
of constant x + is called a front. The FF evolution of a 
quantum system from one front to another is generated 
by the Hamiltonian P~ . 

Despite existence of these options, neutrino oscillations 
have only recently been described using a Hamiltonian 
form of dynamics, and only in the IF Q. The IF Hamilto- 
nian description of Ref. Q is based on a slightly extended 
version of the Gell-Mann and Goldberger formulation of 
scattering theory An extension is needed because 
the formal theory in Ref. Q assumes that the scattering 
region is small in comparison to the volume where the 
incoming particle beams are prepared and outgoing par- 
ticles are detected while the neutrino oscillation experi- 



ments involve scattering regions that extend over so large 
a distance between particle sources and detectors that the 
scattering region in them is actually much greater than 
the entire particle acceleration and detection facilities. 
Thus, the neutrino oscillation experiments belong to the 
class of so-called long-base experiments which involve a 
long distance as an element in the measured observables. 

This article discusses the FF of Hamiltonian dynam- 
ics in application to long-base experiments, following in 
the footsteps of the IF description Q. The main reason 
for introducing the FF in description of long-base exper- 
iments, and particularly for application to the neutrino 
oscillation, is threefold. 

The first reason we wish to distinguish is that the FF 
offers a new relativistic interpretation of the neutrino os- 
cillation in terms of interference of amplitudes mediated 
by virtual states of neutrinos. The new interpretation 
is based on the operational definition of how one traces 
evolution of quantum states that is different from the one 
based on the concept of simultaneity in the IF. Namely, 
one uses the same basic principles that are used in the 
IF but the instants of some selected inertial observer are 
replaced with the fronts. Using the FF of dynamics and 
counting x + instead of time t, one operates with results 
of measurements that are correlated with suitably sent 
waves of light. We shall discuss possible choices for how 
the fronts can be defined and how the physical interpre- 
tation of neutrino oscillations depends on these choices. 

The second reason we wish to distinguish for intro- 
ducing the FF in description of neutrino oscillations is 
that the vacuum problem [4] in relativistic quantum field 
theory is posed in the FF dynamics in a different way 
than in the IF. The difference is a subject of broad inter- 
est and scope and the relevant research on the physics 
of the vacuum is too complex and too voluminous to 
fully quote here [5Hl5j. Consequently, we do not dis- 
cuss the vacuum problem. Instead, we observe that the 
problem of generation of masses of particles in the stan- 
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dard model jig], including the masses of neutrinos that 
are apparent in the neutrino oscillation, can be associ- 
ated with non-trivial properties of the vacuum. Since 
the FF of Hamiltonian dynamics differs from the IF in 
the approach to the physics of the vacuum, and thus may 
also differ in its approach to searches for the dynamical 
origin of the neutrino masses, one is motivated to ask if 
the FF could in principle be used to describe the neu- 
trino oscillation at the current level of its understanding. 
This article provides a positive answer using the effective 
theory that is developed staring from the same effective 
Lagrangian density [l7], EH that was also used in Ref. @ 
in the IF of Hamiltonian dynamics. 

The third reason we stress here is that the FF of Hamil- 
tonian dynamics is distinguished from the IF and PF by 
the fact that 7 out of the 10 Poincare group generators 
do not depend on interactions in the FF, while in the 
IF and PF only 6 generators are free from interactions. 
The 7th kinematical symmetry transformation is a boost 
along z-axis. Formally, the 7th symmetry implies that 
a hadron structure appears the same to all observers re- 
lated by a boost along z-axis. This class of observers 
includes the observer at rest in a laboratory, with re- 
spect to whom a hadron is at rest, and the observer in 
the infinite momentum frame, for whom the same hadron 
moves practically with the speed of light. Therefore, the 
FF of dynamics is considered useful for simultaneous the- 
oretical explanation of hadron structure in the context of 
spectroscopy, where the constituent quark model guides 
phenomenology, and in the context of high-energy scat- 
tering, such as in LHC, where the parton model is used 
to describe structure of hadrons. Therefore, a complete, 
future theory of neutrino oscillations in which the cou- 
pling of massive neutrinos with leptons and quarks via 
massive gauge bosons will be fully understood, is likely 
to require a FF Hamiltonian formulation. 

The article is organized as follows. Section [TTI provides 
a brief overview of a typical long-base scattering sys- 
tem that exhibits neutrino oscillations. We focus on the 
example of an experimental setup resembling T2K [l9j . 
Section IIIII recalls the Gell-Mann-Goldberger scattering 
theory with its extension to the long-base experiments 
and FF of Hamiltonian dynamics. Section IIVI describes 
key elements of the calculation and interpretation of rel- 
evant amplitudes and transition rates. Section [V] con- 
cludes the article. Appendix [5] describes derivation of 
the FF Hamiltonian used in the calculations described in 
Sec. ED 



II. OVERVIEW OF THE SCATTERING 
SYSTEM 

The neutrino oscillation can be observed in various ex- 
periments that differ in terms of energies and momenta 
associated with propagation of neutrinos. When neutri- 
nos interact with hadrons and the energies and momenta 
associated with their propagation greatly exceed masses 



of nucleons, the appropriate dynamical theory should in- 
volve quarks. However, if the energies and momenta as- 
sociated with neutrino propagation are comparable with 
nucleon masses, one may prefer to consider an effective 
theory in which neutrinos couple to whole hadrons, in- 
stead of their constituents. 

The FF of Hamiltonian description of the neutrino os- 
cillation is developed in this article in the context of ex- 
periments resembling the T2K experiment [l9| , in which 
the neutrinos have energies and momenta on the order of 
1 GeV. In this case, it is natural to follow Refs. [l7l [l8j. 
and assume that the appropriate effective Lagrangian 
density for describing interactions between neutrinos and 
protons, neutrons, pions, and muons is 

G f 

d = —j= cosi?c M7 Q (! - 75)^ P7a(l - 9Alh)n 

-t^D^ a (l- l5 )^d a ^ +H.c. (1) 

In this density, there appears the muon- neutrino field v^, 
which is a superposition of three fields of neutrinos v% of 
different masses, to^. Thus, = U^Vi- The same 

interaction Lagrangian density was also used in Ref. Q 
to develop the IF Hamiltonian description of the neutrino 
oscillation. 

In the FF Hamiltonian approach developed here, the 
evolution of a quantum system is traced in the parameter 
x + = x° + x 3 = t + z rather than t. In order to define 
the time t and the space coordinate z, one has to choose 
a frame of reference. We first operationally define a fixed 
frame of reference with coordinates T, X, Y, and Z and 
then use it to describe our different choices of the space- 
time coordinates t, x, y, z in the frames in which the FF 
formalism is developed. 

The long-base experiments are carried out on the 
earth, which rotates around its axis and circles around 
the sun. Therefore, strictly speaking, a frame of refer- 
ence with axes of fixed position with respect to the lab- 
oratories such as Tokay and Kamioka docs not define an 
inertial frame of reference that is suitable for developing 
any approach to neutrino oscillations using the concept 
of inertial observers. However, the corrections due to the 
earth non-uniform motion are small. 

The specific frame of reference with space-time coor- 
dinates T, X, Y, and Z is fixed to the laboratories used 
in a long-base experiment. We shall call its coordinates 
the long-base coordinates. Let us focus attention on the 
example of T2K. In this case, the spatial origin of our 
long-base coordinate system is located in Tokay at the 
outlet of the neutrino source. The time T is defined as 
the time of an observer at rest at this point. The Z-axis 
is directed to the muon detector in Kamioka L ~ 300 km 
away. This means that the Z-axis points about 2° be- 
low the horizontal in Tokay in the direction of Kamioka. 
Let the X-axis be directed south and horizontally, and 
the y-axis nearly vertically up, slightly tilted westward. 
The angular precision of determination of these axes di- 
rections corresponds to the ratio of the size of Kamioka 
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Kamioka Tokay 



FIG. 1: Example of the FF coordinates defined using the long- 
base coordinates as (t,x,y,z) — (T, X,Y, Z). This choice is 
called our preferred choice in the text. The points A and 
B have coordinates x\ = 0, x\ — — 2L, x\ = (0,0), and 
ij = 2L, Xg — 0, xg = (0,0), respectively, where L is the 
long-base distance. The thick vertical lines represent world 
lines of laboratories located in Kamioka and Tokay 



detector (on the order of 30 m) to 2irL, which is about 
0.3'. 

The reference frame (t, x, y, z) of an observer who con- 
structs the FF of Hamiltonian dynamics can be chosen 
in many ways. We particularly distinguish the choice in 
which t = T, x = X, y = Y and z = Z, see Fig. [TJ We 
call this choice our preferred choice of the FF frame since 
it leads to a simple interpretation of the results obtained 
using the FF of dynamics. We shall discuss also other 
choices of (t, x, y, z) with respect to (T, X, Y, Z). 

In the FF with our preferred choice of z = Z, the scat- 
tering system evolves in x + , which is measured along the 
corresponding axis. This axis coincides with the world 
line of a gedanken photon sent at time T = from 
the long-base reference- frame origin in Tokay to Kamioka 
along the if-axis. 

The evolution of the system in x + is generated by 
the FF Hamiltonian, denoted by P~ . Its interaction 
term corresponds to the Lagrangian space-time interac- 
tion density of Eq. ([!]). The Hamiltonian is derived in 
Appendix [Al as an integral of the corresponding Hamil- 
tonian density over the front x + = 0. The FF Hamilto- 
nian density is not as simply related to the Lagrangian 
density as in the IF because of the constraint equations 
that are specific to the FF. In the effective theory we 
use here, the constraint equations, see Eqs. (|A6|) to (|A9|) . 
are coupled and cannot be solved easily. Fortunately, in 
the case of weak interactions, one can use an expansion in 
powers of the coupling constants Gp and . For our dis- 
cussion of the neutrino oscillation, it is sufficient to solve 
the constraint equations using expansion up to terms of 



order GfF^. The resulting Hamiltonian contains all the 
terms that are required for a description of the leading 
effect of neutrino oscillation and provide its FF interpre- 
tation. The generic interpretation involves the following 
experimental setup. 

The 7r + -beam is assumed to move nearly along the 
Z-axis, a bit upward in the long-base coordinate sys- 
tem. Let the p, produced in decays of 7r + move at even 
greater angle upward. In such cases, the physical four- 
momentum transferred from Tokay to Kamioka via in- 
termediate quantum states, p v = (E v ,p£ ,p% ,p„), may 
in the majority of muon detection events in Kamioka be 
reconstructed as lying close to the four- vector with com- 
ponents [E v , 0, Q,p%). This means that in our preferred 
frame of reference the physically likely four-momentum 
transfers have components close to the four-vector p u — 
(E u ,0,0,p*) with pi — pZ. The corresponding FF mo- 
mentum four- vector coordinates are p^j — E v ± p* — 
E v ± an d Pu — 0- Of course, the FF components of 
p v would be different if the front x + = were chosen in 
a different way. For example, if one reversed the z-axis, 
the FF components would be p^ — E v ^p^ and pj; = 0. 

The FF formal scattering theory describes the prepa- 
ration of the 7r + beam in terms of its gradual buildup in 
x + rather than t. The beam-preparation time r intro- 
duced in Ref. Q following Gell-Mann and Goldberger Q 
in the IF of dynamics, is replaced by the parameter r + 
in the FF of dynamics. Using the convention that c = 1 
and e" = 2/r+, the factor e T ' T = e eT in Refs. @, Q is 
replaced below in the FF of formal scattering theory by 
the factor e x+/r+ = x+ / 2 . Note that X+ denotes the 
auxiliary FF "time" variable over which one integrates 
in analogy to the integration over the time variable T in 
Eq. (2.3) in Ref. @. 

Physical interpretation of the pion beam preparation 
in the FF of Hamiltonian dynamics depends of how one 
chooses the z-axis in the long-base frame of reference. In 
our preferred FF frame, the buildup of the initial pion 
and neutron states has a relatively simple interpretation. 
The pions emerge from a carbon target under "pressure" 
exerted on it by highly energetic protons. The target 
has a shape of a rod, its transverse cross-section diame- 
ter being much smaller than its length. The uncertainty 
of created pion position (momentum) is comparable with 
the size (inverse of the size) of the carbon target, denoted 
by r (from the word "rod" ) . With the use of electromag- 
nets, pions produced at various positions with various en- 
ergies and moving at various angels are assumed focused 
to travel approximately in one direction. This increases 
the intensity of the pion beam and results in some energy 
distribution of pions in the beam. One could attempt to 
describe this distribution using a carefully adjusted den- 
sity matrix. In this paper, we limit our consideration 
to the case of a mono-energetic pion beam with a well- 
defined momentum. 

We can estimate an upper limit on the value of a FF 
beam-preparation "time" r + by considering the pion life- 
time. Namely, the length of the tunnel in which 7r + moves 
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after being produced and in which it is nearly certainly 
turned into /2 and by the electroweak interactions is 
on the order of 100 m. Since the pions travel with nearly 
speed of light, they cover approximately the same dis- 
tance in time and in z-direction in our preferred frame. 
Thus, the estimate for the upper limit on r + can be taken 
as 2 times the length of the tunnel. This length is much 
shorter than the long-base length L. Thus, the approx- 
imate, formal Gell-Mann and Goldberger description of 
the pion wave function buildup ought to contain a fac- 
tor e x+ / T+ , where r + <C L. Then, e~ = 2/r + is much 
greater than 1/L. The value of e~ that corresponds to 
at most 100 m is on the order of at least 10 -9 eV. At 
the same time, the expected differences between the FF 
neutrino "energies," p~ , with masses squared on the or- 
der of 10~ ( 3 ^ 7 ) eV 2 divided by p+ ~ 1 GcV, are on the 
order of 1(H 12 - 16 ) eV, i.e., they are much smaller than 
e~. This means that in our preferred frame the FF "en- 
ergy" eigenvalue density of final eigenstates of P~ has a 
width much larger than the differences among eigenval- 
ues of Pq~ of virtual neutrinos with different masses rm in 
the intermediate states. This fact will be shown to lead 
to the validity of the standard oscillation formula in the 
FF of Hamiltonian dynamics. 

If we chose the z-axis in the opposite direction, i.e., 
z = —Z instead of z — Z, the carbon target, the pion 
tunnel, and the proton and pion beams would all ap- 
pear nearly instantaneous in x + = t + z. However, the 
strength of the pion wave function would not immediately 
be normalized to 1 in the quantization volume according 
to the Gell-Mann-Goldberger formulation of scattering 
theory. Namely, following the Gell-Mann and Goldberger 
model for the beam-buildup process in the IF, one can 
postulate that in the corresponding FF model the wave 
function strength should smoothly increase with x + no 
matter how the z-axis is chosen. In the case of z = —Z, 
one considers the incoming plane- wave state |</>i) whose 
probability increases with X + according to the function 
g £ x+/2 f or x+ < q f or vames f x - anc [ x ±- "simul- 
taneously" in the sense of x + . 

Now, when z = — Z, the duration of the whole scatter- 
ing process in the sense of x + is also short. This happens 
despite that the corresponding "distance" x~ between 
the pion source and final muon detector is very large, 
about twice the distance L between Tokay and Kamioka. 
In fact, it will be shown below that x + ~ {p+ / E V )L, 
where p+ is very small for z = —Z . In Sec. lIVl we will dis- 
cuss the corresponding quantum-mechanical mechanism 
by which the same oscillation formula works irrespective 
of how the FF z-axis is introduced. 

In the FF, we shall encounter instantaneous interaction 
terms, in the sense of no delay in x + . These interaction 
terms are called seagulls. They correspond to the so- 
called Z-diagrams in the IF of dynamics in the infinite 
momentum frame (IMF), originating in virtual particles 
with momenta oriented opposite to the infinite momen- 
tum. Since the meaning of "simultaneity" in the sense 
of x + depends on how x + is defined, which in turn de- 



pends on the choice of the z-axis, we will discuss the role 
of seagulls in neutrino oscillations for different choices of 
z-axis. The same physical phenomenon of oscillation will 
result from considerably different accounts of the flow 
of x + and correspondingly different accounts of the in- 
termediate quantum states that depend on the different 
choices of z-axis. 



III. FRONT FORM OF SCATTERING THEORY 

The IF version of scattering theory needed for de- 
scription of long-base experiments has been obtained in 
Ref. [2| using the Gell-Mann and Goldberger formal ap- 
proach as a starting point Q. The same principles are 
followed here to obtain the FF version of required theory. 
In the case of T2K, the theory objective is to calculate the 
transition rate from the state \<f>i) = \tt + ti) to the state 
\4>f) = \pfJ-p) as a function of distance between the pion 
source and muon detector (effects of the Fermi motion are 
ignored here, see comments in Ref. [1] and Sec. E]). The 
formalism does not depend in its essence on the states 
details and it is developed using general notation of \(j>i) 
and \<pf) for initial and final states, respectively. 

According to the discussion in Sec. [Ill the incoming 
state is gradually built up according to the formula 

|*^ + )> = f dX+e^ x+ ' 2 

' J — oo 

xr ,p-(x+-^)/2| $i ^ +)) (2) 



The state \§i{X + )) is a function of X + that results from 
action of exp (— iPq X + /2) on \<j>i). The operator P ~ is 
the FF free Hamiltonian, such as in Eq. (|A14I) in the 
case of T2K, and the state \<pi) is its eigenstate, with the 
corresponding eigenvalue denoted by p~ . 

The transition rate of the evolving system to the final 
state \<j)f) can be measured in terms of counting parti- 
cles in detectors (such as Kamiokandc) using x + or using 
t. It is natural for physicists to use t. This is the time 
of an observer at rest in our long-base reference frame 
(Kamiokande is at rest with respect to this frame). Us- 
age of x + appears less natural. However, the FF count- 
ing of transition rates is physically quite realistic in the 
sense of counting particles in coincidence with detection 
of light that defines the fronts of varying x + . In the long- 
base experiments, where the size of a neutrino detector 
is negligible, one can take advantage of the fact that the 
entire detector world line has a fixed value of z. Thus, the 
measurement of a FF transition rate in terms of states lo- 
calized in the detector corresponds to the differentiation 
of detection probability with respect to x + and 

d dt d dz d 

~dx~+ ~ 'dx+dt + dx+'d'z' ^ ' 
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For detectors located at fixed positions in the long-base 
frame of reference, the operational definition of transi- 
tion rates in x + is obtained by setting dz/dx + = and 
observing that dt/dx + = 1/2. 

The probability that the system is in state \<j>f) at x + , 

is 



UJfi(x + ) 



\A(x+)f 



(5) 



where the scattering amplitude calculated following the 
Gell-Mann-Goldbcrgcr formulation of scattering theory 
is 



le e 



A{x + ) = (0/1 — 



i(p~ -p-)x+ /2 



pT - p - 



(6) 



and the norms of states are not changing with x + . In 
the ratio of transition rates we are going to calculate a la 
Ref. @, , the norms cancel out and they are omitted in 
further discussion (an alternative is to think about them 
as equal 1). Using identity Q for differentiation of uifi 
with respect to t, and omitting the norms, one obtains 



d_ 



\A(x + )\ 2 



dx+/2 



\A(x+)f 



(7) 



for detectors located at fixed z. Following the steps anal- 
ogous to the IF calculation @|, one obtains the FF ex- 
pression for the transition rate, 



dx+/2 



\A(x + )\ 2 



2e 



(P f -P] ) 2 + (e~) 



■\Rfr{x + )\*, 



(8) 



where PL (x + ) is 



Pi — P + ic 



Mi 



(9) 



This result is used in the next section to derive the neu- 
trino oscillation formula. 



IV. FRONT FORM OSCILLATION FORMULA 

There is a difference between the IF and FF calcu- 
lations of transition rates that results from a difference 
between the Hamiltonians obtained from the same La- 
grangian density of Eq. (fT]) in these two forms of dynam- 
ics. The difference between the Hamiltonians is a conse- 
quence of the constraints that appear in the FF and are 
absent in the IF of dynamics. The FF interaction Hamil- 
tonian density is not merely a negative of the Lagrangian 
interaction density (with removed time derivatives of the 
pion field). Namely, it only depends on the dynamically 
independent components of fermion fields and contains 
terms that are instantaneous in x + , called seagulls. 



Appendix [S] describes the calculation of P~ . The 
constraints in effective four-fermion theories, such as in 
Eq. (JTJ) , are a set of coupled non- linear equations. We 
can only solve them using a perturbative expansion in 
powers of Gf and F n . Fortunately, including terms or- 
der Gf, F w , and GfF^ is already sufficient for describing 
the dominant effect of neutrino oscillations. 



Calculation of R 



fi 



Denoting terms order Gf and as Pf , and terms 
of second order including terms order GfF^ as P 2 one 
obtains from Eq. (J9j) that 



Pi" 



O i(p--P )x+/2 



Pi 



Pn 



-Pl+P 2 



(10) 

\4>i), 



where the phase factor in front is not important because 
of the modulus in Eq. ([8]). The first term in the square 
bracket describes transitions with a neutrino or an an- 
tineutrino in the intermediate state. The second term 
comes from the seagull interaction. 

We denote by p u the four-momentum — pp that is 
physically transferred from 7r + to n. The transfer is car- 
ried by neutrino or anti-neutrino or mediated by the seag- 
ull term. The four-momentum transfer does not depend 
on the kind of neutrino that appears in the intermedi- 
ate state. The seagull term sums up effects coming from 
constraints on neutrinos of all masses. 

A priori, there are two possible types of the interme- 
diate states, both coming in a sum over neutrino kinds 
labeled using subscript i. For p+ > 0, the intermediate 
states contain p,, n, and i/£. For < 0, the interme- 
diate states contain n + , p, fx, and Pi. No matter which 
possibility one considers for incoming and outgoing par- 
ticles, the ones for which the reconstructed p + > or the 
ones for which p+ < 0, there is always an additional con- 
tribution from a seagull. The sign of reconstructed p+ 
may a priori also depend on the choice of z-axis. How- 
ever, when the experimentalists focus on events in which 
the reconstructed momentum transfer p v is obtained as- 
suming that p^ > 0, only states with virtual neutrinos 
can contribute and anti-neutrinos are excluded. This 
is the case when the leading contributions to the total 
scattering amplitude come from the reconstructed four- 
momentum transfers p v that must be close to an on-mass- 
shell four-momentum of a neutrino with some non-zero 
mass. These are the cases we focus on here. We shall 
come back to the possibility of studying anti-neutrino 
oscillation effects in Sec. IIVDI 

In our preferred frame, every neutrino carries a large 
positive p+ . The phase factor in the exchange term reads 



exp{i( Pi -P )x + /2} = exp{i{p u -p )x + /2} 



(11) 
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Using the interaction Hamiltonian of Appendix \K\ the 
matrix elements that occur in Eq. (|10[) are obtained in 



the forms 



(pmmI-Pi — 



Pi - p o + ie 



PT\mr + ) = -igfSipt ~ P 1)5 {2) {pi ~ p})u pl a {l - g Al *)u n 



pt Pu - Pvi + it 



(pfifi\p- f \nn + ) = -igfS(pt- P j)S^( P i-pj)u pl a (l-9Al 5 )ur l 



x u^ la (\ - 7 5 )-l T (_ l ^) 7/3 (l _ 7 5 ) v 



(12) 



(13) 



r 



where P g j obtained from the density of Eq. (| A16I) is the 
only term in P 2 ~~ that contributes. 

B. Neutrino exchange vs. seagull 

Our goal now is to show that the oscillating exchange 
term in Eq. ([10[) dominates in the transition rate. 

The two contributions in the square bracket in Eq. 
(1T0|) . shown in Eqs. ([12]) and dHJ), differ only by the fol- 
lowing factors that appear between the two current fac- 
tors that are expressed in terms of the spinor matrix ele- 
ments associated with the pion-decay and neutron-proton 
transition: 



1 



pi Pv 



Pvi 



2 P i 



(14) 



These two terms can be looked at as contributing 7- 
matrices with various coefficients. The seagull con- 
tributes only to the coefficient of j + . One can see the 
relative size of the two terms by first factoring out the 
Feynman-like denominator, 



1 



1 



1 



Di{p v ) pi(p v - Pui +ie ) pl-m^.+ie pi' 

(15) 

in both of them, and then comparing (assuming pi = 0) 



Pu 



(Pu -P Vi + i€ pi")-^--(!6) 



The second term, that comes from the seagull, is as small 
as the off-shellness of the neutrino, and the first term 
contains pi. and p~. of which at least one is large; which 
one depends on the choice of z-axis. The FF smallness 
of the seagull term in comparison with the neutrino ex- 
change term corresponds to the IF smallness of the anti- 
neutrino exchange terms in comparison with the neutrino 
exchange terms in the Hamiltonian description of neu- 
trino oscillations in Ref. [2J. 



Thus the exchange term with an oscillating phase fac- 
tor is expected to dominate the transition rate. We have 
inspected matrix elements of the above matrices in ran- 
domly selected momentum and spin configurations that 
may be relevant in experiments such as T2K. In all cases 
we so sampled, our inspection confirmed that the above 
described analysis of the coefficients of 7-matrices cor- 
rectly estimates the actual ratio of the full amplitudes 
and the seagull contributions can be neglected in com- 
parison with the exchange term contributions. Once 
the seagull term is neglected, the sum over intermediate 
states with different neutrinos with appropriate phases is 
inserted in Eq. ©. 



C. FF interpretation of the oscillation 

In a theory of the quantum mechanical scattering pro- 
cess, one can fix momenta of the initial and final parti- 
cles. These momenta define the momentum transfer four- 
vector p v . The physical transfer of p~ can greatly differ 
from p~. , since the latter is calculated from the neutrino- 
mass-shell condition for given pi and pi , while the for- 
mer results from the difference of momenta of scattering 
particles. However, due to the denominator Di(jp v ), the 
scattering amplitude is large only for p~ near one of the 
values p~. . In experiments such as T2K, it may be as- 
sumed that the observed counts of muons in the far detec- 
tor come from these large terms. Since the actual size of 
the scattering amplitude also depends on the uncertainty 
of the FF "energy" e~ , the final interference pattern be- 
tween the amplitudes coming from specific intermediate 
states also depends on the size of e~. 

Fig. [2] shows that the three amplitudes that contribute 
to counting muons coming from interactions that involve 
three different intermediate states with different-mass 
virtual neutrinos can interfere in analogy with the in- 
terference pattern that is familiar from the elementary 
quantum slit interference experiment, except that the 
role of slit position is played by the FF on-mass-shell 



"energy" p~ . The number of "energy slits" available for 
the interference depends on the size of c~. which must be 
large enough to enable the interference of all potentially 
available states with virtual neutrinos. 




Kamioka 



Tokay 



FIG. 2: Schematic representation of the interference of am- 
plitudes mediated by virtual states with neutrinos of different 
FF free "energies" pZ i , due to the different masses m^. The 
Pu ( increase with the mass m,; p vi < p V2 < p v;i . Like in the 
elementary quantum slit experiment, amplitudes with certain 
values of p~ provide dominant contributions to the interfer- 
ence pattern. The FF "energy" uncertainty e~ must be large 
enough to create the standard interference pattern in the total 
counting rate of muons in the far detector. 

For x + = 2L, the neutrino mass-dependent phase fac- 
tors present in the complete scattering matrix element 
Rji {x + ) in Eq. (|8|), can be reduced to simpler forms, 
exp(ip^.2L/2) — > exp(im^.L/p+), since the transverse 
and + momenta are the same for all kinds of the con- 
tributing neutrinos and the common phase factor due to 
p^ is irrelevant to the interference pattern. Thus, the 
amplitude R% (x + ) in Eq. ([Tllj) significantly varies with 
p v only due to the identified above mass dependent phase 
factors and FF "energy" denominators for different neu- 
trinos. The result for R^ (2L) is proportional to 



E 



u, 



in 



Di{p v ) 



E 



U, 



(17) 



The statement made above that the interference of ampli- 
tudes with different neutrinos in the intermediate states 
occurs for sufficiently large value of the uncertainty e~, 
is illustrated in Fig. [3] When the size of e~ exceeds the 
differences among the on-mass-shell FF "energies" p~ , 
all neutrino "channels" or "slits" in p~ contribute and 
the corresponding total amplitude exhibits the standard 
interference pattern. 

On the other hand, a reduction of the ratio of the ex- 
perimental p~ uncertainty, e - , to the FF "energy" dif- 
ferences among different neutrinos, would lead to devia- 
tions from the standard oscillation formula, since some 
"channels" or "slits" in p~ would no longer be able to 
contribute to the lepton counting in the far detector. For 
example, heavy sterile neutrinos with sufficiently large 
masses must drop out form the interference pattern with 
relatively small uncertainty e~. 



a) 






FIG. 3: Qualitative plots of the function | 2t=i V-^*(P")I us ~ 
ing arbitrary units for a) e~ = 0, and b) e~ 7^ 0. The neutrino 
masses m Vi are assumed to increase with i; m vi < m„ 2 < 
m„ 3 . The scales are arbitrarily adjusted for the purpose of 
obtaining a clear drawing that explains how it happens that 
the larger e~ in comparison to A) 
the standard oscillation formula; 
the more neutrino intermediate states coherently contribute 
to the muon counting rate in the far detector. 



iijL/p^ the more accurate 
the wider the width e~pt 



It is clear that the experimental and theoretical ways 
of studying deviations from the standard oscillation for- 
mula require further investigation. Namely, the patterns 
of deviation depend on the neutrino masses (e.g., see Ap- 
pendix in Ref. !2j) and coupling constants, including the 
factors Ufj,i. Therefore, the deviations are sources of in- 
formation on these parameters. The way to reduce the 
relevant ratio or, equivalently, to enhance the FF "en- 
ergy" differences in comparison to e~ using the front 
choice of z = Z, is to reduce pj. In order to have a 
possibility of reducing p+ and still creating the final lep- 
ton, one might consider detection of electrons instead of 
muons. 



D. 



Other choices of z-axis 



In the experimental T2K-like setups where the ob- 
served neutrino interference patterns are unambiguously 
identified using reconstructed four-momentum transfers 
p v , one can establish if the physical p v must approxi- 
mately match the value that corresponds to a 7r + -decay 
into some free neutrino and p,. If it does, Vi and p, may 
only form states with a free invariant mass comparable 
to the 7r + mass and with a relatively small uncertainty 
determined by the product e~pj. In such cases, the p„ 
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is always greater than no matter how one chooses the 
z-axis to define the front. All intermediate states of the 
virtual particles that can significantly contribute contain 
only neutrinos. The seagull contribution is always small. 

Still, the theoretical FF interpretation of the oscilla- 
tion changes when one changes z-axis. A particularly 
instructive example of a change is provided by the case 
of z = —Z, see Sec. [TTJ In this case, assuming p^j ~ 0, 
the FF "time" interval x + during which the neutrino ex- 
change occurs, estimated as the time difference between 
the duration of a flight of a real massive neutrino and a 
massless photon from Tokay to Kamioka, approximately 
equals (p+ / E V )L. For a neutrino on the mass shell m,-, 
one obtains x + ~ xf ~ (mf /2E^)L, instead of 2L ob- 
tained in the case z = Z discussed in the previous sec- 
tion. But the corresponding neutrino's "energy" is now 
approximately p~ = 2E U . Therefore, the half of a prod- 
uct of a short x\ and a large p~ that appears as a phase 
in the exponent is the same as the result mfL/(2E [/ ) in 
the case z = Z . The width of Feynman-like denomina- 
tors is also obtained without change. Namely, although 
the p^, for z = — Z is very small, the pion beam is pre- 
pared nearly instantaneously in x + in comparison to 2L. 
This means that e~ is large and the product e^pj is not 
changed. In summary, the intermediate neutrinos that 
are nearly instantaneously exchanged in x + for z = —Z 
have very small p^ but a huge p~ and hence the oscilla- 
tion pattern is not expected to change by changing the 
front despite that the FF interpretation depends on the 
choice of a front. The interpretation certainly differs from 
the commonly applied imaginative picture of a neutrino 
that flies from place to place changing its mass "state." 

Besides the two choices of z = Z and z = —Z, there 
exists a whole set of choices with z-axis directed at some 
angles to the Z-axis. In principle, one can also discuss 
reference frames in motion with respect to the long-base 
frame. Instead of discussing these multiple options, we 
prefer to point out a qualitatively different aspect of the 
FF interpretation of the neutrino oscillation. 

Namely, one may ask if the four-momentum transfer p v 
that is reconstructed in analyzing observed events must 
always be close to a p v that corresponds to a real decay 
of 7T+ or, perhaps, it could involve p+ < 0. In the latter 
case, the cc + -dependent phase factor for z = —Z could be 
interpreted as due to a free evolution of the anti-neutrino 
that nearly instantaneously in x + travels from Kamioka 
to Tokay, where it is absorbed by a pion. The authors 
have not succeeded in finding out if even in principle the 
existing or prospect data on the neutrino oscillation per- 
mit a reconstruction of p v without assuming that it must 
approximately correspond to a real pion decay, instead 
of absorption of an anti-neutrino that nearly instanta- 
neously in x + is created in Kamioka and subsequently 
absorbed in Tokay. It is also not clear if such reconstruc- 
tion, even if possible, could be associated with a mea- 
surable counting rate of muons in Kamioka. The only 
observation that we can offer is that the invariant mass 
difference between a neutron and a proton-muon pair im- 



plies that the reconstructed p v would have to correspond 
to an anti-neutrino virtuality on the order of 100 MeV, 
which suggests greatly reduced muon counting rates. 

On the other hand, the intriguing aspect of consid- 
ering the reconstructed p+ close to is that the FF 
vacuum problem can be associated with the region of 
vt = [H [HI [H|- Thus, the opportunity for recon- 
structing momentum transfers with extremely small p+ 
in long-base neutrino oscillation experiments could per- 
haps be used to test theoretical ideas concerning vacuum 
involvement in the generation of neutrino masses. Stud- 
ies of the small p+ region would have to include a precise 
description of the Fermi motion effects associated with 
the binding of a neutron in the oxygen nucleus in the 
far detector (in the case of water detectors) and with the 
binding of quarks inside hadrons. In any case, the region 
of p+ = is singular in the quantum field theories that 
form the standard model and deserves a formal study in 
the context of neutrino oscillations. 



V. CONCLUSION 

The FF of Hamiltonian dynamics provides an alterna- 
tive description of the neutrino oscillation to the formal 
description available in the IF of Hamiltonian dynam- 
ics [2| . Since the FF also provides an alternative formula- 
tion of the ground-state problem in quantum field theory, 
including the concepts of vacuum condensates and mass 
generation, we conclude that the neutrino oscillation can 
be studied using the FF of dynamics in conjunction with 
the fundamental issues of particle theory. 

The FF approach provides an interpretation of the 
standard neutrino oscillation formula as resulting from 
the interference pattern that occurs only when the ex- 
perimental "energy" uncertainty is sufficiently large 
in comparison to the differences between the individual 
on-mass-shell values of p~ for any of the neutrinos of the 
kind i. The FF explanation of the interference leads to 
the condition that when the differences between individ- 
ual p~ , for the same p+ and p-^ , are much greater than 
the uncertainty e~ , the standard oscillation formula is 
not valid. While the conditions of validity of the standard 
oscillation formula are well satisfied in the case of experi- 
ments like T2K and 3 already known neutrinos with quite 
small masses, they are not satisfied for much heavier neu- 
trinos. 
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Appendix A: Calculation of P 

The operator P~ is a generator of translations in x + . 
It is defined in the FF of dynamics through an integral 
over x~ and x , 



p- = i / ,/.,•- ,/-.,- J T 



(Al) 



where T^ v is the energy-momentum density tensor and 
half of T+- will be denoted by V~ . 

Standard methods of canonical quantization [22| lead 
to the expression for P~ that corresponds to the La- 
grangian density C = Co + £j, where 

£o = d^d^n - ml^n + ^2$(iB- m^ijj ,(A2) 



and £7 is given in Eq. ([T]), according to the formula 

TtW = 9" a Edd^ d ^-^ C - ( A3 ) 

<t> 

Summation over <fi is meant to indicate that one sums 
over all fields in the theory. is expressed in terms 

of the fields and their conjugated momenta in order to 
quantize the theory. 

The canonical conjugate momenta for the fermion 
fields, d£/dd~ip = ipij + /2, depend on the field compo- 
nents ip^> = ^7°7 + V': which are the dynamically inde- 
pendent variables. Matrices A^ = ^7°7 ± are projectors. 
The fermion field components = ^7°7 _ V' satisfy 

constraint equations which couple all fields in the theory 
through the interactions. Constraint equations for 
follow from the Euler-Lagrange equations, 



(A4) 



dd a ip dip 

which can be written for every fermion field ip in the form 

{i$-m^)il> = id+j ^, (A5) 

where ipi denotes the interaction terms. Namely, 
1 



Vil 



: U; t T a ^(gnT a A p-zfd a ^) , (A6) 



id+ 

/</ --^T^ (gpT"n + if d a 7r) , (A7) 



id 
1 



Pi 



(A8) 
(A9) 



where = 7 Q (1- 3A7 5 ) and P = 7 «(l- 7 5 ). One can 
define a free field 



(A10) 



where ipQ ' = ip( ) — ip\ ' and ip^ = ip( +s> satisfy 
dition 



con- 



4>V = + 4 +) ■ (All) 



A calculation yields 
dC 



dd~ip 



d ip — ip(i$ — m^)ip 



(A12) 



^07 



+ -(d ± ) 2 +m 
2id+ 



-^ -iJiA-(td + )ip I 



The interaction parts of all fields can be expanded 
in a series of powers of the coupling constants g — 
G F cos6» c /\/2 and / = The resulting FF Hamil- 

tonian density V~ takes the form of a series 

V~ = Vq + V{ + V2 + o(g k f:k + l<2), (A13) 

where 



</>0 



2id+ 



-V>o 



(A14) 

= (i//2)3VJ M - ifd^Jto - gJ^Juaa + H.c. , 

(A15) 



and denotes all terms order g , f and gf. The 
subscript in currents and J£ indicates that the nu- 
cleon current Jjy = P7 Q (1 — gAl$)n and lepton current 
J£ = J2i ^ui^»7 Q! (l — 75 )m are evaluated with all the free 
fermion fields that are generically defined in Eq. (|A10|) . 

It turns out that among all terms in V% only the 
terms proportional to gf are important for calculation 
in Sec. IIV1 These are the seagulls, 



■P-j = -*5//Wot(1-7 5 )^M1-7 5 )W) 



21 
2id+ 
+ 



-igfv^{l - 7 5 )^+ /jvo(1 - 7 5 Ko 



+H.c. 



(A16) 
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